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Abstract 

The cubic anisotropy model provides a simple example of a system with an 
arbitrarily weak first-order phase transition. We present an analysis of this 
model using e-expansion techniques with results up to next-to-next-to-leading 
order in e. Specifically, we compute the relative discontinuity of various phys- 
ical quantities across the transition in the limit that the transition becomes 
arbitrarily weakly first-order. This provides a useful test-bed for the applica- 
tion of the e expansion in weakly first-order transitions. 
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I. INTRODUCTION 



The cubic anisotropy model is a simple two-scalar model that, for a certain range of pa- 
rameters, has a phase transition with similarities to the finite-temperature phase transition 
of electroweak theory in the early universe. More generally, it provides a simple proto- 
type for systems that have weak, fluctuation- induced, first-order phase transitions In 
ref. [^], we and our collaborators discuss in detail the similarities and dissimilarities with 
the electroweak transition (noted earlier by Alford and March-Russel |^) and described how 
the cubic anisotropy model is a good testing ground for analytic techniques that claim to 
distinguish between second-order and weakly first-order phase transitions. In particular, the 
model can be used to study e expansion methods, which we have previously applied to the 
electroweak case f^. This paper presents the details of calculating e expansions for weakly 
first-order phase transitions in the cubic anisotropy model. A better overview of the moti- 
vation, a summary of our e expansion results, and a comparison against numerical Monte 
Carlo simulations may be found in ref. 0. 

Our goal will be to compute the relative discontinuity of various quantities (the specific 
heat, susceptibility, and correlation length) across the phase transition when one has an 
extremely weak first-order transition. Specifically, we will compute ratios such as x+/X- 
where x± ^i-re the susceptibilities on either side of the transition. This was originally done at 
leading order in e by Rudnick some twenty years ago We have extended the calculation 
to next-to-leading order (NLO) for the correlation length and next-to-next-to-leading order 
(NNLO) for the susceptibility ratio. A next-to-leading order calculation for the specific heat 
ratio is performed in a companion paper 0. 



Our leading-order results for x+lx- ^-^d however, differ by factors of 4 from ref. 
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A. Three-dimensional reduction and the cubic anisotropy model 



Before introducing the cubic anisotropy model, we shall very briefly review the connection 
between phase transitions in thermal quantum field theory and those in classical statistical 
mechanics. For definiteness, consider the topical example of electroweak theory. In studying 
the electroweak transition, one starts with a 3+1 dimensional SU(2) gauge-Higgs theory at 
finite temperature.^ Schematically, the Euclidean action is of the form 

S = I'dr I d'x!^lm^ + \F^-^\M' + + ■■■], (1.1) 

(with gauge-fixing terms omitted). (3 is the inverse temperature, and c? = 3 is the number 
of spatial dimensions. If the correlation length at the transition is large compared to the 
inverse temperature (which is generally the case), one may simplify the study of equilibrium 
properties of the transition by integrating out the dynamics of the Euclidean time direction. 
This yields an effective three-dimensional theory that describes the long distance physics of 
the transition and which may be precisely matched, order by order in coupling constants, 
to the original theory:^ 

Se, = p J d'^x!^lm' + \F^ + l\<p\' + ^\M' +■■■]■ (1.2) 

For a review, see refs. [§],§[. There is no need to go into detail here, except to note that the 
mass-squared t of the Higgs in the effective theory has the form 



^Fermions, and the U(l) and SU(3) gauge fields, do not have a major impact on the phase transition 
dynamics and are, for simphcity, neglected. 

^ Another way of explaining the appearance of a three-dimensional theory is to note that, if effective particle 
masses are small compared to T at the transition, then for small momenta the Bose distribution function 
l/(e'^'^— 1) is large compared to one. But physics should be classical if the number of quanta in each state 
is large. The long-distance physics of the transition can therefore be approximated by classical statistical 
mechanics in three spatial dimensions. By the well-known equivalence of statistical mechanics and quantum 
mechanics, this is equivalent to a "zero-temperature" field theory in three Euclidean space-time dimensions, 



which is one way to view the effective theory (1.2). 
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t = -fi' + cg'T' + -- - , (1.3) 

for some constant c. The fact that t becomes positive as the temperature T increases drives 
the restoration of manifest SU(2) symmetry at high temperature. The action ( p..2[ ) describes 
a classical statistical mechanics problem in three spatial dimensions, where Ses is to be 
interpreted as PH. 

In this work, we will not study the three-dimensional theory ( |1.2| ), but will instead 
examine a simpler three-dimensional theory consisting of two scalar fields known as the 
cubic anisotropy model. In a more general form, the cubic anisotropy model is an 0(n) 
symmetric scalar model of n real scalar fields, to which is added an interaction that breaks 
0{n) symmetry down to hyper-cubic symmetry IHJ^. The action is 



S = I d'^x U 



2 t 

+ 2 



+ KU 



:i.4) 



(Note that the overall (3 in ( p..2| ) can be absorbed by a rescaling of 0.) We will ultimately be 
interested in the simplest case, n = 2. The parameters u and v are dimensionless coupling 
constants, and k is a dimensionful normalization which we will fix later. The phase transition 
of interest occurs as the parameter t is varied. At tree level,^ the transition appears to be 
second-order, with hyper-cubic symmetry spontaneously broken for t < and restored for 
t > 0. As we shall review, however, the effect of higher-order corrections on the nature of 
the transition cannot be ignored. 



^ There is some ambiguity of language depending on whether one views the action ( |l.4| ) as (a) describing 
classical statistical mechanics of a field theory in d=3 spatial dimensions, with S equaling f3H, or (b) as 
a quantum-field theory in d—3 Euclidean space-time dimensions. In the former case, a "tree level" result 
would normally be referred to as a "mean field theory" result; in the latter, it would be referred to as a 
"classical" result. The first interpretation more accurately reflects the physics of the problem, but the latter 
is more familiar to particle theorists. We shall bypass the terminology issue by referring simply to tree-level 
vs. one-loop results, etc. 
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B. The e expansion 



In second-order phase transitions, the correlation length diverges at the transition and 
physics near the transition is dominated by large infrared fluctuations which cannot be 
treated perturbatively. One way of summarizing this for a theory like the cubic anisotropy 
model is to note that, since the couplings ku and nv of the theory have non-trivial mass 
dimension A—d = 1, the dimensionless loop expansion parameter R must, by dimensional 
analysis, be of the form 



which diverges with the correlation length. In first-order phase transitions, the correlation 
length is finite, but perturbation theory can still fail if the correlation length is large enough 
that R ^ 1. We shall refer to such a situation as a weakly first-order transition. In this 
paper, our goal will be to study the arbitrarily weak limit R oo. 

The e expansion is based on generalizing d=3 spatial dimensions to 4=4— e dimensions. 
When e is small, one can systematically remedy the problems of perturbation theory by 
using suitable renormalization-group (RG) improved perturbation theory. Computing to 
some order in RG-improved perturbation theory corresponds to computing to some order in 
e. At the end of the day, one sets e— i>l in the resulting truncated series. In some cases this 
is known to give quite good results.0 



^ A slightly more detailed review for particle theorists, in the context of the electroweak phase transition, 
may be found in the introduction of ref. jsj. 



R 




X 
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FIG. 1. Renormalization group flow into the infrared of {u,v) for small e and (a) n < rtc = 4 + 0(e), 



or (b) n > ric- The lightly shaded region is the domain in which the tree-level potential (1.4) is unbounded 
below. The more heavily shaded "wedges" show the regions in which perturbation theory is reliable. In the 
left-hand figure (a), the special trajectories labeled A and B will be a focus of attention. These are the limiting 
trajectories for theories which approach the Ising or cubic fixed points arbitrarily closely before flowing off 
toward the classical instability line. Such theories have arbitrarily weak first order phase transitions. 

Fig. shows the renormahzation group flow, for small e, of the dimensionless couplings 
u and V of the cubic anisotropy model as one moves to longer distance scales. The lightly 
shaded region, delimited by the lines u = —nv and u = —v, designates the range of couplings 
where the tree-level potential of ( |1.4|) is unbounded below. For t> = 0, the cubic anisotropy 
model reduces to an 0{n) model, which has a second-order phase transition associated with 
the infrared fixed point marked "0(n)". For u = 0, the model ( |1.4| ) reduces to n uncoupled 
copies of a basic quartic scalar field theory, which is in the same universality class as the 
Ising model. The associated fixed point is marked "Ising" in fig. |I[ Besides the Gaussian 
fixed point at u=v=0, there is another fixed point known as the cubic fixed point. All of 
these fixed points occur at couplings of 0(e), and so the couplings may be treated as small 
and perturbative when e ^ 1. 

The stability or instability of these fixed points depends on the number n of scalar fields. 
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For n < He = 4: + 0(e) (which encompasses our main case of interest, n=2), the 0(n) fixed 
point is infrared stable,^ as shown in fig. |l|(a), while the Ising and cubic fixed points have 
an unstable direction and correspond to tricritical points.Q A theory whose couplings lie 
between the lines running from the origin to the Ising and cubic fixed points, respectively, 
will fiow at large distances to the 0(n) fixed point and so will have a second-order transition 
with 0(n) symmetric critical behavior. A theory with couplings outside of this region 
(and not on the critical lines bounding this region) does not fiow to any weakly coupled 
infrared-stable fixed point, and so might be expected to have a first-order phase transition. 
This is indeed the case. As discussed by Rudnick [|l|, and as we shall review below, the 
criteria R <^ 1 for the success of perturbation theory (known as the Ginsburg criteria) is 
satisfied for couplings very close to the region of tree-level instability, designated by the 
heavy shaded regions in Fig. |l]. In these R <^ 1 regions, one indeed finds that perturbation 
theory reliably predicts a first-order transition. The renormalization group fiow then shows 
that any theory whose couplings are outside the boundary of the basin of attraction of the 
0{n) fixed point is equivalent to a theory with couplings having i? <^ 1, and so will have a 
first-order transition. The lines from the origin through the Ising and cubic tricritical points 
are therefore boundaries separating theories with first- and second-order transitions. 

The case of n > ric is shown in fig. |l|(b). The 0(n) and cubic fixed points exchange roles 
relative to the n < ric case. 

Now consider a sequence of theories with first-order transitions, such as those indicated 
by the dotted line near the top of fig. |l|a, which approach the Ising line m=0. The correlation 
length at the transition can then be made arbitrarily large, since for m=0 it is infinite. The 
renormalization group trajectory approaches the dashed line in fig. ^ labeled trajectory A. 
This limiting trajectory first fiows into the Ising fixed point along the line m=0, and then 

^ Except, of course, with respect to the parameter t (i.e., temperature), which is not shown in the figure 
and which has to be fine-tuned to reach the phase transition. 

Further expansion in e yields = 4 - 2e + 2.588e^ + O(e^) §]. 
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flows away from tlie Ising fixed point in the unstable direction toward the region of classical 
instability. A similar limit for the tricritical behavior of the cubic fixed point gives trajectory 
B. This is the limit we will take to obtain arbitrarily weak first-order transitions in the cubic 
anisotropy model. For n = 2, trajectories A and B are equivalent because a redefinition of 
by a 45° internal rotation, 

{(f)i,(f)2) + (f)2,<f)i- h) , (1-6) 
leaves the Lagrangian ( p.. 41) in the same form but with 

iu,v) ^ iu + ^v,-v) . (1.7) 
This means that theories below the u axis in fig. |l|(a) are, for n=2, related by the mapping 



LTD to theories above the axis. 

We will therefore compute properties of the transition, for small e, on trajectory B. This 
is easiest to do by following the trajectory into the perturbative region R <^ 1. Throughout 
this paper, we focus on the flow away from the cubic fixed point. For n=2, the final results 
for physical quantities must be the same as for flow from the Ising flxed point. For other n, 
the flow from the Ising flxed point could be analyzed similarly, but we have not bothered to 
do so. 

In the next section, we flx notations and renormalization scheme conventions. In sec- 
tion |IT|, we review the leading-order analysis of the susceptibility ratio x+/X-i which was 
originally carried out by Rudnick . The most straightforward derivation uses a calculation 
of the one-loop effective potential and the explicit one-loop RG equations. After reviewing 
this calculation, we show that explicit knowledge of the one-loop potential and RG equa- 
tions was not actually necessary. In section |^ we extend the calculation to next-to-leading 
order. This calculation does require the explicit one-loop potential and RG equations, but 
we show that explicit knowledge of the two-loop corrections (which would be used in the 
most straightforward derivation) is not required. Section |V] computes the next-to-leading 
order result for the correlation length ratio ^^/^^. Then, in section |V|, we flnally ex- 
tend our calculation of x+ 1 X- next-to-next-to-leading order, which requires a non-trivial 



ring-diagram resummation of perturbation theory, the exphcit two-loop potential, and the 
explicit two-loop RG equations. Finally, section |V11| summarizes our results. A review of 
the leading-order result for the flow of the couplings {u,v), originally derived by Rudnick 
as well as details of the two- loop potential, are left to appendices. 

II. NOTATION AND CONVENTIONS 

We will use dimensional regularization for loop calculations in d = 4— e dimensions and a 
renormalization scheme closely related to modified minimal subtraction (MS). Specifically, 
the bare Lagrangian is 

/:bare = |^||50r + V^bare(0) (2.1) 

with^ 



Zuu\(l)\^ + Z^v^i 



(2.2) 



and 

Ar={47c'f^r[l-i), (2.3) 



and where all renormalization constants have the form 



e2 



(2.4) 



Note that we have rescaled our couplings and A by an additional factor of (47r)^ compared 
to the typical convention in particle theory. The additive constant A is irrelevant to the 
calculation of the susceptibility or correlation length ratios and may be ignored, but it will 
be important for the specific heat ratio computed in ref . . Note that we have relabeled the 
parameter t as m^, which is the typical notation used in particle theory. But it should be 
kept in mind that variation of really represents variation of temperature in the underlying 
physical problems of interest. 



^For n — 2, our couplings (u,v) are related to Rudnick's 0| choice of couplings, call them (uR,fR), by 
UR — {u + v)/l2 and wr, = ~v/6. 
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The susceptibility x is defined by adding a linear term hcj) to the Lagrangian and defining 

X-lim^t*. (2.5) 

h^o dh 

The normalization and regularization of the hip term is not important because we will ulti- 
mately only be interested in the ratio x+lx-i where it cancels out. 



III. REVIEW OF LEADING-ORDER ANALYSIS 

Two things are needed to analyze the transition: (a) the location of some point (m, v) 
along the portion of the trajectory within the perturbative regime, and (b) a perturbative 
analysis at that point. Technically, it is easiest to choose the point where the trajectory 
intersects the line u = —nv or u = —v, where the tree-level potential first becomes unstable. 
(The full, effective potential remains stable, as it must, since it is invariant under changes 
of renormalization scale.) We shall call this point {u^:,v^,). 

The value of {u^:,v^:) turns out not to affect the leading-order calculation of x+/X-i so 
we shall proceed for the moment without it. (The value of the couplings does affect the 
leading-order results of other ratios, such as the specific heat ratio C+/C_ IHJ^I-) 

At tree level the transition appears second-order. For > 0, the minimum of the tree- 
level potential is at 0=0. For < 0, it is along an "edge", oc (1, 0, 0, ■ ■ ■), if f < 0, and 
along a "diagonal", oc (1, 1, 1, ■ ■ if > 0. Rather than discussing the entire structure 
of the effective potential V{(j)), it will generally be sufficient simply to consider its behavior 
in the relevant direction. For flow from the cubic flxed point {v < 0), attention can be 
restricted to an edge. The tree-level potential Vq then becomes 

K)(0) ^ K)(0) = ^^202 + ^^^^^^{u + v)^' , (3.1) 

where 

0= (0,O,O,---). (3.2) 

At the instability line, u=—v and hence the quartic interaction term disappears along the 
edge. 
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To see the first-order nature of tlie transition, one must consider the effect of the first 



loop correction. This is just the Coleman- Weinberg effect |T0|. For the sake of definiteness, 
and because we need the resuhs later on, we shall go through the explicit calculation of the 
one- loop corrections to the effective potential. However, after the fact, we shall show that 
an explicit calculation was actually unnecessary for computing at leading order. 



A. The one-loop potential 



The one-loop contribution to the effective potential is 



Vi{(f)) = X(m^) + (ra— l)X(m^) + (counter-terms) 



(3.3) 



where and are the eigenvalues of the curvature of the tree- level potential Vo(^) 
evaluated at 0, and the one-loop integral I{z) is 



I{z) = |trln(-92 + z) = -\ {A7t)-^/^t[-I) 



A/2 



(3.4) 



Along an edge, we have 



ml = m^ + lJ\ffi^{u + v) 4^- 



(3.5a) 
(3.5b) 



which are the curvatures parallel and orthogonal to the edge, respectively. 

Now take u = —v to fix ourselves on the tree-level instability line. It is notationally 
convenient to express the potential in terms of 



M2 = i/xW' 



U I 



(3.6) 



and one finds 



XlJ.'Vo{(l)) = A + 3u-^m^M^ , 



m 



2e c/(rf-2)sin(7re/2) 



Ye rf(rf-2)sin(7re/2) 



(3.7) 
. (3.8) 
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For the moment, we're only working to leading-order in e, so we can take the limit e— >0 to 
find:R 



Afli'{Vo + Vi)=A + Su^^m'^M^ + \^^ 



In 



In I 



m 



+ 0(e). 



(3.9) 



As one varies m? , this potential describes a first-order transition which occurs at an 7^0.0 
If it weren't for the one-loop corrections, the transition would be second-order and occur at 
m^=0. Since couplings are 0(e), one then expects that at the transition is small if e is 
small. We shall indeed see a posteriori that 



— ~0(e) 



(3.10) 



when the order parameter is the same order of magnitude as its value in the asymmetric 
phase. So in this range of we can drop m compared to M and find 



M^l'{Vo + VI) = A + ?,u-^m^M^ + 
where 



+ 0{tM^,m'M\m^), (3.11) 



C 



n—1 
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C 



10 



Cii 



(3.12) 



The above approximation to the potential has two degenerate minima when 



2 2 

m = rrii 



Cii u jj? 



■ exp 



C 
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3 " \ Cn 
and the value of M at the asymmetric minima is 



n—1 
~V1 



(3.13) 



^ Note that leading order in this context doesn't simply mean the tree-level potential. We are interested in 
the leading-order results for quantities describing the first-order nature of the transition. But the tree-level 
potential by itself does not describe a first-order transition. 

"'^"By transition, we mean the point where the two ground states are degenerate. In physical applications, 
this may not be the point of direct physical relevance if there is significant super-cooling. 



13 



M2 = Ml = /i^exp f-1 -^)= _ (3^^4) 

V Cii / 

As promised, rn^/Mf ~ 0{u) ~ 0{t). The curvatures at the origin and at the asymmetric 
minima are0 



= |Ar/iW'(0) =m2[l + 0(e)], (3.15) 

0=0 

= iAr/i^Ml/"(Mi) = 2m2 [1 + 0(e)] (3.16) 



(where V\M) = dV{M)/dM, etc.). Thus 

X+ 



--2 + 0{e). (3.17) 

X- 



B. Avoiding one-loop details 

In the preceding derivation, the final result for x+/x~- did not, in fact, depend on the 
values (|3.12| ) of the constants Cu and Ciq. So we could have arrived at the same result 



simply knowing the form (|3.11|) of the one-loop potential in the limits u=—v, and 

■m? /M"^ ~ 0{t). But this form just follows from (a) the existence of the renormalization 
group, which produces a single power of In /i at one-loop order, and (b) the fact that M is 
the only relevant dimensionful parameter if m is negligible. The logarithm must therefore 
be ln(M/ /i) in our approximation. If we had understood a priori the independence of 
the result on Cu and Cio, we could have avoided doing the explicit one-loop calculation. 
Generalizations of the following arguments will later save us from the need for two-loop 
calculations at next-to-leading order, and three-loop calculations at next-to-next-to-leading 
order. 



The proportionality relationship reflects the fact that the derivatives on the right-hand side are with 
respect to 0(/i) — the renormalized field at the scale where w(/i) = —v{fi). However, this (p is proportional 
to the bare 0, and thus the ratio (a) does not depend on the proportionality constant, and (b) is insensitive 
to short-distance physics. 
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The independence from Ciq may be understood as follows. The only parameters that 
final results can depend on are the dimensionless couplings {u^:,v^) and the corresponding 
renormalization scale fi. Other parameters, such as and the scalar expectation M have 
been solved for and eliminated by requiring that we be at the transition and in one or the 
other phase. If our final result is a dimensionless ratio, such as x+/X- or C+/C^, it must 
then be independent of /i and can depend only on So the answer can't change even 

if we arbitrarily change fi in (|3.11|) while holding ('U*,f*) fixed. (This is different from a 



simple statement of RG invariance, which would involve changing {u, v) in a compensating 
manner when changing /i.) Such a change in ^ (at lowest order) is equivalent to varying 
Cio, and so dimensionless ratios do not depend on Ciq. 

The coefficient Cu of the logarithm is determined by the one-loop renormalization group. 
However, as noted above, it is not needed for x+/X-- The reason is that, at leading order, 
the result (|3.17|) is independent of the couplings {u^:,v^). So the result would be the same 
if we redefined (u, v) by (m, v) {xu, xv) for some constant x. This redefinition does not 
take us off the line u=—v, but it does change the coefficients of the one-loop /5-functions and 
therefore changes Cu. The leading-order result for x+/X- must therefore be independent 
of Cii. This simplification doesn't occur for C^/C_, which turns out to be proportional to 

at leading order 

As a prelude to our later higher-order analysis, it will be useful to sketch the renormal- 
ization group determination of Cu. The one- loop RG equation for the effective potential in 
the cubic anisotropy model is 

where^^ 

l3uiu,v,e) = -eu + /3u{u,v) , (3.19a) 



"'^^The trivial e dependence in the beta functions (3.19) is, of course, a standard feature of minimal sub- 
traction renormalization schemes. 
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f3v{u,v,e) = -ev + (3y{u,v) , (3.19b) 

with 

^u{u, v) = f3^J^ + 0{u\ v"^) , = u (|(n+8) u + 2v) , (3.20a) 

^^{u, v) = + 0{u^ v^) , =v{4u + 3v) , (3.20b) 

P^2{u, v) = Pill + 0{u^ v^) , = |(n+2) u + v, (3.20c) 

-f^iu,v)=Oiu^,v^). (3.20d) 

For the susceptibihty, we are only interested in the (j) dependence of the potential, and the 
running of A will be irrelevant. As it turns out, we will also not need (3^'^ or 7^ below. 

The renormalization group flow does not map the line u+v = onto inself. Consequently, 
to apply the RG equation, one must retain the v dependence of the tree-level potential ( p.l| ) 
rather than specializing from the outset to u=—v. Working at e=0, one easily finds that the 
RG equation is satisfied by 

V, + Vr = V,-{\ /3« m V' + ^ Ar[/5(^) + /^(i)] 0^ } In /i 

+ (/i-independent) + 0(eV) . (3.21) 



Taking u=—v (and neglecting m? relative to M^) then gives 



(3.22) 



which agrees with ( |3.12| ). In later sections, this same RG method will be used to determine 



//-dependent terms in the two- and three-loop contributions to the effective potential. 



C. Scale hierarchies and subtleties at higher orders 

As we shall see, the preceding tricks for simplifying calculations will generalize to higher 
orders in e as well. Leading-order results for generic ratios require only explicit knowledge 
of the tree-level potential and the one-loop renormalization group, and x+/X- doesn't even 
need the latter. Next-to-leading order calculations generically require the explicit one-loop 
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potential and the two-loop renormalization group, although x+/x~ needs only the one- 
loop renormalization group. Next-to-next-to-leading order generically requires the two-loop 
potential and the three-loop renormalization group, and so forth. 




FIG. 2. A sequence of diagrams of the same order in e in the asymmetric phase. Sohd lines represent 
heavy (wb) degrees of freedom; Hght hnes represent the light (ma) degrees of freedom. 

There is a subtlety, however, in the simplistic assumption that each successive order 
in e requires exactly one more order in the loop expansion of the effective potential. The 
source of this subtlety is the ratio of scales m? /M"^ ~ 0{e) in the asymmetric phase. Fig. |^ 
shows a sequence of diagrams, with arbitrarily many loops (starting at 3 loops), which are 
all the same order in e in the asymmetric phase. These diagrams consist of multiple "ring" 
corrections to the light (mass m?) mode due to interactions with the heavy (mass M^) 
modes. The outer loop is dominated by momenta of order m, and the cost of adding an 
additional ring is 

O('uV'M') = 0{uM^/m') = 0(1) . (3.23) 

This particular problem could be handled diagrammatically by resumming the light 
propagator to incorporate all one-loop heavy rings. A more elegant way to think about it is 
in the language of the renormalization group. At distances large compared to 1/M in the 
asymmetric phase, our scalar theory should be replaced by an effective theory consisting of 
only the light degree of freedom. The mass of the light scalar in the effective theory will be 
its original mass m plus corrections from integrating out the heavy modes. 

The free energy in this light effective theory will be of order m^. Comparison to ( ^.11[ ) 
then reveals that it will only be important at next-to-next-to-leading order. If only working 
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to NLO, one can ignore the need for this resummation in the asymmetric phase. 



IV. NLO ANALYSIS OF x+/x 

A. The effective potential: asymmetric phase 

Going one order beyond the previous analysis requires consideration of (a) two-loop 
contributions to the effective potential, and (b) corrections to the e — and m? /M"^ — > 
limits we took of the one-loop potential. For the latter, we can simply expand the general 
one-loop potential (|3^ ) to the desired order: 



M^i'{VQ + Vl) = ^ + 3u~^m^M^ + 
with 



+ Ar/iWi + 0(eV)asym, (4.1) 



2Cnlnf— j -2Cn 



elO 



and 



a 



elO 



(4.2) 



(4.3) 



The subscript "asym" is a reminder of the assumption tv? /M"^ ~ 0(e) in the error estimate, 
which is valid only in the asymmetric phase. 

For the two-loop contribution, we may take e=0 and ignore m? altogether. The renor- 
malization group requires the contribution to have the form 



M2 



M2 



C.ln^ - +C.ln 



G 
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+ 0(eV) 



asym • 



(4.4) 



As in the leading-order calculation, we will not actually need to compute all three parame- 
ters. 

The first simplification is to note that rescaling /i by 



H ^ nil + xu) 



(4.5) 
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in Vq + Vi + V2, while holding all couplings fixed, changes C20 at this order but nothing else. 
Therefore dimensionless ratios cannot depend on C20 at NLO. Similarly, a change such as 



/X (1 + xe) 



(4.6) 



for a NLO 



would change C^io and nothing else, so we never actually needed its value 
calculation. 

We could determine all the other constants in ( |4.4| ) by requiring the potential to satisfy 
the RG equation at two loops. However, analogous to what happened at leading order, we 
will not need all of these coefficients for x+/X-- It is sufficient to apply the RG equation 
at one loop order, as given by ( |3.18| - pl20D . However, we do need the one-loop potential for 
general {u,v) without the restriction u=—v. Returning to ( p.3| ) and ( p.5| ), one easily finds 



AA/i^(K) + K) = A+- 



3 u+v ^4 , 9 {u+v f 



2 



4 

I (n-1) 



In 



In 



u /i^ 

'M2\ _ 3 



3 

) 2 
+ 0(eV) 



asym 



(4.7) 



Applying the one-loop RG equation]^ and then setting u=—v, determines 

n+2 



c. 



22 



6 



-C 



11 • 



(4.8) 



It is also worth noting that it was unnecessary to compute explicitly the 0(e) correction to 
Vi given by the first term of ( [4. 21) , because the coefficients of the logs in that correction are 
determined by the RG equation as well, arising from the explicit e in (3.19). This observation 
will substantially simplify the analogous calculation when we later proceed to NNL order. 

To find the asymmetric phase susceptibility X-^ we now need at next-to-leading order 
both the value of at the transition and of in the asymmetric phase. Perturb around 
the one-loop solutions ( |3.13|) and ( |3.14| ) by writing 



V = 1^(1) + 6V 



m 



ml + 6m'^ 



M = Mi + 5M, 



(4.9) 



"'^^When applying the renormalization group equation, it is helpful to note that is multiplicatively 
renormalized, ii{dM^/dfi) = AP Pm^, with = Pu/u + 27^ = + 0{u^). 
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where V(i) is the one-loop e=0 approximation (|3.11 



) to the effective potential and 



5V = 6Vi + V2 



(4.10) 



as parameterized in ( |4.2| ) and ( [4.4[ ). By linearizing the equations V {M)=V (0) and V'{M)=Q 
that determine the asymmetric phase expectation M at the transition point, one finds 



and 



6m 



6M 



Ml, mi 



M2 d (W 



Ml, mi 



The fractional shift in 1/x- is 
Putting in the explicit form for 5V , we obtain 



' — ^ + — + ^TTTT 



Ml ,mi 



m = m 



= Ml 



13 C, 



elO 



16 C 



(3n+2) 3 C 
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20 



l/X- oc 9^1^(0) = 2m^ 



11 
1 

1 - -e 
2 



2Cii Cii 



(n-13) 



12 



u + Ofe' 



(4.11) 
(4.12) 

(4.13) 

(4.14) 
(4.15) 
(4.16) 



Note how all the dependence of d'^V on C^io, C21, and C20 in the result ( [4.16| ) is hidden in 
the overall factor of m?. 



B. The effective potential: symmetric phase 



When examining the asymmetric phase, we made an expansion in m^/M^ = 0(e). For 
the symmetric phase, where M=0, this is not a good approximation. The symmetric phase 
is easier, however, because one needs only one-loop contributions at e=0 for next-to-leading 
order results. Consider the one-loop potential (|3.9| ). It gives 



1/X+ oc aJ\/(0) = 1 + 



<^ 1 



(n-1) 
6 

(72-1) 

6 



u 



u 



In 



m 



\ 12 , 



(4.17) 
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Putting this together with the asymmetric phase result gives our NLO ratio 

X+ 



X- 



. 1 {n-l) , Mn-l) \ 
1 e in + 



(4.18) 



where we now explicitly remind ourselves that u is to be evaluated at the point {u^,v^). At 
leading order, we didn't need to know at all for x+/X-- For next-to- leading order, we 
need to find the leading order value of m*. 

C. {u.t,v.f) at leading order 



The one-loop RG equations for the couplings are 



(4.19a) 
(4.19b) 



where (5^'^ and are given in ( p.20| ). The explicit solution was found when n=2 by 
Rudnick [|1| and generalized to other n by Domany, et al. We review the derivation in 
Appendix ^. The resulting trajectories are given by 



where 



i?(/,c) = 
A = 



V = eR{u/v, c) , 



{nX + l){nX + 2) 
3 



(nX + 1 



X / -2/ + - 
A n 



-nX 



4 - n 



(4.20) 

(4.21) 
(4.22) 



and each choice of the constant c picks out a different trajectory. The trajectory that fiows 
away from the cubic fixed point is c = 0; the one fiowing away from the Ising fixed point is 
c = 1. The values of {u,v) on the tree-level instability line are then 



, , , , r./ N S{n^ + 5n + 3) 

(m*, v^) = (m,, -M*) , = -eR{-l, 0) = e , 

n[n+2){n+6) 



(4.23) 



for fiow from the cubic fixed point and 
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3n 



(n+2)(n+8) 

for flow from the Ising fixed point. For n=2, these are 



3- (1(1-^; 



(7n+2) 



e, (4.24) 



VSO"' 80^) '^'■^^ \ 160"' SO^J ' 



(4-25) 

respectively, and are related by the mapping ( [1.71) . 

Our final, next-to- leading order result for x+/X- is given by (|4.18| ) and ( [4.231) . For n=2 
this yields 



X+ 
X- 





a7e\ 


\l eln( 


L 160 ^ 


.3207 



(4.26) 



The presence of an e In e term is a feature that one does not encounter in e expansions for 
critical exponents of second-order transitions. It arises here from the hierarchy of scales 
m^/M^ ~ e characterizing the physics of the asymmetric phase at the transition. 



V. NLO ANALYSIS OF ^l/^l 

The correlation length ^ is determined by the location of the pole of the two-point 
correlation. This is the solution p"^ = —^^^ to 



+ + n(p2) = , 



(5.1) 



where n(p^) is the (one-particle irreducible) self-energy. Since m^+n(0) is another name for 



the susceptibility 1/x, we can write 

= + fn(-r') - n(o) 



(5.2) 



As we shall see, this equation can be solved by iteration, treating the second term, call it 
An, as small. Hence, at leading order is the same as ^ind is 0{m'^). Fig. |^ shows 
the only one-loop graph that contributes to the momentum dependence of n(p^). 



22 



FIG. 3. One loop diagram contributing to momentum-dependence of 11 in asymmetric phase. 



In the symmetric phase, the expectation (0) is zero and so fig. ^ vanishes. The only mass 



scale is m and so two-loop contributions to AH would be order 0{v?'m?) = 0(e^^~^), which 
can be ignored at next-to-leading order. 

In the asymmetric phase, the largest correlation length will be that associated with the 
degree of freedom 0i along the edge, corresponding to = O(m^) in (|3.5|) . For u = —v, 
this degree of freedom does not couple to itself, and the degrees of freedom running around 
the loop in fig. |^ are the heavier ones (pi {i ^ 1) associated with the mass scale ml = 0{M'^). 
The momentum dependence AH of fig. ^ is therefore 0{up^) = 0{uim}) = 0(e^~^). As 
advertised, it can be treated as a perturbation. Explicit calculation in the -C limit 
gives 



n(p') - n(o) 
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asym ) 



(5.3) 



and so 




(5.4) 



(5.5) 



Combining with the result ( [4.18| ) for the susceptibility ratio gives 




■u, +0(e'). 



(5.6) 



Finally, inserting the n=2 value (^4.25|) of u^, yields 




(5.7) 



23 



VI. NNLO ANALYSIS OF 



A. The 2-loop renormalization group 



Two-loop RG /^-functions for the cubic anisotropy model may be easily extracted by 
following standard derivations in one-scalar models and replacing the overall couplings of 
each diagram by those appropriate for our two-scalar model. One finds 



with 



(2) 



(2) 



Pm2 



(2) 



(3n+ 14) , 22 2 5 2 

u U V uv , 

3 3 3' 



46 



17 



(5n + 82) 2 2 - 3 

u V uv V 

9 3 3 



(n+2) 2 9 



{n+2) 
36 



-u 



-uv 
6 



12 



(6.1) 
(6.2) 

(6.3) 
(6.4) 
(6.5) 
(6.6) 



B. The effective potentiah asymmetric phase 




+ 




FIG. 4. One-loop diagram contributing to effective light mass. 



As discussed in section [III C| , the NNLO analysis of the effective potential in the asym- 
metric phase requires separating the light and heavy modes and doing a resummation of 
the heavy modes' effect on the light ones. The effective mass nies of the light mode at dis- 
tances large compared to 1/M can be found by explicitly computing the diagrams of fig. ^. 
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More simply, it can be taken from the curvature of the one-loop potential near the 



asymmetric minima ( |3.14D : 



m;s{M) =m^ + 3- -M^ + (n-1) uM 



u 



'J MI 



(6.7) 



2 V /i- 

The sub-leading 0{eni}) corrections to the above relationship are convention dependent: 



they depend on exactly how we want to define meg. However, such sub- leading corrections 
to m^g are not relevant at the order of interest. 

The effective potential in the asymmetric phase is 



light 



^hcavy 



r2T/heavy\ 



asym ; 



(6.8) 



where Vq + Vi^^''^ and 5Vi'^^^^ are given by ( [4.1| ) and (^4.2|). We have added the super- 
script "heavy" to indicate that these contributions come from the heavy modes. The next 
correction ^^^^heavy ^YiQ expansion of the second term in the one- loop potential 



IS 



/ A/f2 



12- 



emlO 



■^Cnln^g) +|C,oln^(^) -laioln(^) + C.^o 



. (6.9) 



This can be justified either by explicit expansion or by application of the renormalization 
group. The new constant C^eio will be irrelevant because it can be absorbed into a redefinition 
of /i. The new constant Cemio may be extracted from explicit expansion of (|3.9| ): 



a 



emlO 



1 + ^ 



(6.10) 



The last piece of the one-loop potential we need is the contribution from the light modes, 
corresponding to the third term of ( |3.9| ). However, this contribution must be computed with 
the correct effective mass ( |6.7| ) so that 



ET/light 



1 



eff 



In 



m 



eff 



3 
2 



(6.11) 
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As usual, this mass resummation is most easily accomplished by rewriting the light mass 
term of the Lagrangian as 



treating the first term on the right-hand side as part of the unperturbed Lagrangian and 
the second term as a perturbation.]^ This perturbation will generate a new graph at NLO, 
shown in fig. ^ 



FIG. 5. Next-to-leading order contribution involving the effective mass counterterm (m'^—m'^g). 

For the two- loop potential, first consider the e — > and m —>■ limit of (|4.4| ), which we 
now refer to as V2'^'^'^^. For the current calculation, we will need to know all of the coefficients 
{C2i }. C21 may be determined by either explicit calculation or by applying the two- loop 
RG to the one- loop potential (^?7|). The constant C20, however, requires explicit calculation. 
The two-loop contributions are given the diagrams of fig. combined with fig. ^ 




(6.12) 





X 



FIG. 6. Two-loop diagrams contributing to V2 



The details of the calculation are given in Appendix |^, and one finds 



^21 — 7: — L'li , — tv: 

6 6 



(6.13) 



See, for example, ref. |12|. 
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in addition to the previous result ( ^l8|) for €22- 

The sub-leading corrections to V2^^^^ come from relaxing the e^O limit and remembering 
that the heavy mass ( |3.5bD is + instead of simply M^. Expanding in e and m? /M"^ 
gives 



Ciiln^ 



^ ^ C'lllnf — + Cm20 



(n+2) 
12 



Cnln^^l 



'ri+10 



3 

(3n+14) 
12 



+ 



e20 



(6.14) 



V 3 36 / V /i^ 

where the coefficients of all the logs can be determined by requiring the full potential to be 
invariant under the two-loop RG. The constant C^2o "will be irrelevant as it can be absorbed 
into a redefinition of /i. The remaining constant, 6*^205 is found by expansion of the explicit 
(e=0) two-loop potential, as described in Appendix One finds 

1 



a 



m20 



c 



11 



(6.15) 



Finally, the contribution of heavy modes to the three-loop potential has the form 



C33 InM ^ ) + C32 In^ I ::- 1 + C31 In ( 1 + C30 



+ 0(eV)asym. 

(6.16) 



We will only need C33 and C32 for x+/X-- (C'si can be changed by a suitable redefinition of 
the couplings, and so cannot affect the physical ratio x+lx--) The coefficients C33 and C32 
can be determined by applying the two-loop RG to the full potential. To do so, we first need 
to relax the restriction u=—v in our calculation of the two-loop potential. The analysis can 
be simplified a bit, however, in that we can treat u+v as small and only keep terms linear 
in u+v. (That is because the RG derivative (3udu + (3vdv acting on u+v does not give zero 
when u = —v, although the same operation on [u+vY does yield zero.) Keeping only such 
terms in the m=0, approximation to the potential gives the following analog to 



Xfx'iVo+Vi+V2)=A + M' 



C 



10 



+ u 



C^22ln'(^| + C2iln(::- 1 + ^20 
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3 (u+v) . . 



3Ciiln 



+ 0{{u+vf) + 0{e)^. 
Combining with V3, applying the RG, and then taking u=—v yields 

C33 = 3^ (n' + 8) Cu , C,2 = {2n' + 45n + 7) Cn 



5Cn 
2 

(6.17) 



(6.18) 



C. The effective potential: symmetric phase 

In order to find the critical value of m^, we must equate the free energy of the two phases. 
The asymmetric phase approximation ( |6.8D is not good in the M=0 symmetric phase because 
it relies on the approximation m^/M^ ~ 0{e). We must compute the symmetric-phase free 
energy independently. At the order of interest, it is just the one-loop contribution ( |3.9| ): 



n 



J\ffi'V{0) = A + - m^ 



(6.19) 



For the NNLO susceptibility in the symmetric phase, we need the curvature of the 
two-loop potential at M=0. This is obtained by differentiating the general result ( pT]) of 
Appendix ||,[] and yields (lOll) . 



D. {u^,v^) at NLO 

The NNLO result for x+/X- will require a NLO value for (m*, f*). Consider the one-loop 
result for the RG trajectories, ( |4.2(]| ) with c = or 1, flowing from the tricritical points. 
Consider the solution for u as a function of / = u/v, and call it uW(/). Now look at the 
perturbation as we include higher loops. Start by making the rescaling 



fi = fi 



(6.20) 



Replace J and / by J and / in (Bl) and remove 1/e poles, as discussed in the appendix. 



28 



which makes the e expansion exphcit in the RG equations: 

fi— = -V + (3i^\u,v) + e (3i^\u,v) + ■ ■ ■ , 
where /^^"^ is the n-th order contribution to the /3-f unction. Now expand 

uif)^umif) + e5if) + 0ie'). 



(6.21a) 
(6.21b) 



(6.22) 



Plugging into the renormahzation group equations, hnearizing in the perturbation 5, and 
solving yields 



5(/) = e^(^) y^^d/'e-^(^') 



+ e^(^)5(/o), 



(6.23) 



(/')>/' 



where 



K(f) = f df 



d 




du 




u 




^2 



"111 (/'),/' 



and we have defined 



(6.24) 
(6.25) 

(6.26) 



For the trajectory flowing from the Ising flxed point, we should find that the dependence on 
the initial perturbation 5(/o) vanishes when /o approaches the Ising tricritical line /o=0. In 
this case, the system will first flow to the Ising flxed point before flowing away, thus washing 
away dependence on S{fo)- Similarly, for flow from the cubic flxed point, the dependence on 
S{fo) should vanish as /o — > —A. This is indeed the case. For any trajectory, one flnds 



(6.27) 



which vanishes as /o ^ or —A. 

The resulting correction 5{f) does not seem to have a simple form for general n, but 
may be evaluated explicitly when n = 2. Taking /o = —A and / = — 1 for the cubic point 
trajectory, one flnds 
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(6.28) 



One may verify that the analogous calculation for the Ising fixed point trajectory yields the 
appropriate transformation (|1.7| ) of ( |6.28|) . 



E. Results 

By equating the potential in the symmetric and asymmetric phases, one determines 
and the asymmetric phase value of to next-to-next-to-leading order, on the line u+v = 0. 
One finds 



= ml < 1 + e 



and 
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+ eu 



47n 247 7(n+6) 2 C'.2o 
+ -rr H tt:: — 



384 64 



288 



+ 



1152 96 
+ 0(e^)+0(u')] 



144 



6 



144 
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2 Cii Cii 

(6.30) 



The curvatures of the potential in the two phases, at the transition, when u+v = 0, are 



dlV{0) =mUl + u 



(n-l) 
6 

n-l) 



In 



n—1 
12 ' 2 



e u 



(n-l) 
24 



-In^ 



n—1 



n—1 



-u + In u + 1 

12 / V 12 ' 



36 



n+2) In 



9 / n—1 \ , / n—1 
2' M -91n 



12 



12 



u 



12v/3(^ln2-.(0)-3 

(6.31) 
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dlV{<P) = 2rr? 



e — n ■ 



(n-13) 



2 ~ " 12 

(n-l)(27-2n) 



eu ■ 



{n-9) 
12 



72 

+ 0(e3) + 0(M=^)|. 



12 



72 



(6.32) 



Here, L{z) is Lobachevskiy's function, defined in Appendix H and 



L( ^ ) = 0.02461715- ■■ . 



(6.33) 



As required, all dependence on the undetermined parameters C30, C31, etc., is hidden in the 
overall factor of m^. The susceptibilities x+ ^-^^d X- equal these curvatures up to a common 
overall proportionality constant. Inserting the value of derived in the previous section 
yields our final result for the susceptibility ratio at n=2, 



^ = 2 1 
X- I 



iIlnflIi)+H 
160 V320/ 80 
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1602 



In^ 



/17e\ /8374 
V320/ ^ VT6O2 
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VII. DISCUSSION 



We now collect our results for n=2 and evaluate the coefficients numerically: 



^ = 2 

X- 



^2 

2± — 9 



1 + e(-0. 1063 lne + 0.4494) 
+ e2(-0.00731n2e + 0.1647 In e - 0.2859) + 0{e^] 
1 + e(-0. 1063 lne + 0.4139) + 0(e^ " 



(7.1) 



(7.2) 



The ratio C^/C^ of specific heats will be evaluated in ref. 0, and all of these results are 
compared against Monte Carlo simulations in ref. 0. 



31 



Looking solely at the series above, our results are moderately encouraging. At e=l, 
the NLO corrections are 45% for x+/X- ^"^^ 41% for NNLO corrections for x+/X- 

drop to 29%. The sub leading corrections could have turned out quite large compared to 
the leading-order results, as is known to happen in some cases where the leading-order e 
expansion is a poor quantitative approximation. Q In contrast, our series seem tolerably well 
behaved. 

It should not be too difficult a calculation to extend ^+/^?_ to NNLO, but we have not 
done so. An interesting question for further research is whether it is possible to determine 
the large-order behavior of e expansions for the ratios we have investigated. The techniques 
used for critical exponents of second-order transitions]^ do not obviously generalize to 
this problem. 

This work was supported by the U.S. Department of Energy grants DE-FG06-91ER40614 
and DE-FG03-96ER40956. We thank David Broadhurst and Joseph Rudnick for useful 
conversations. 



APPENDIX A: REVIEW OF LEADING-ORDER {u^,v^ 



Consider the one-loop RG equation given by (|4.19| ) and ( p.20|) . Before we look for 



trajectories, note the location of the fixed points at leading order in e: 



Gaussian: 


(u. 


v) 


= (0,0), 


Ising: 


(u. 


v) 


= (o,t), 


OH: 


(m. 


v) 


~ ( (n+8) ' ^) ' 


cubic: 


(m. 


v) 





See, for example, the discussion of large n in ref. 

For a review, see sees. 27.3 and 40 of ref. and references therein. 



(Al) 
(A2) 
(A3) 
(A4) 
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Though we shall not directly make use of it now, note that the dependence on dimension in 
the one-loop RG equations can be eliminated by the rescaling ( |6.20| ) so that ( [4.19[ ) becomes 



(A5a) 
(A5b) 



This depends on the fact that the (3^^^ are quadratic in u and 



The following is a sketch of the solution to the RG equations (|4.19|) as found by Rudnick 



HI] and Domany, et al. [|Tl|]. Begin by removing the first term on the right-hand side by 
switching to new variables A and 5, 



so that 



u = jjL "A^ V = ji . 



(A6) 



(A7) 



Divide these two equations, and note that the right-hand side is a function solely of / = 
u/v = A/B: 

dA (3W(A,B) 



H{f). 



(A8) 



dB f3i'\A,B) 

Changing variables from {A,B) to {f,B), solving the resulting equation, and fixing the 
boundary condition -B(/o) = -Bo gives 

rf df 



B = Boexp (^J^ 
With the beta functions ( |3.20| ), we have 



(A9) 



/i u 
fi" V 
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(fo + \\ 















nX 



(AlOa) 
(AlOb) 



with A given by ( [4.22|) . Now, from f = u/v, note that 
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Of! V 



(All) 



Use the solutions (|A10|) to write the right-hand side in terms of /, /i, and constants. Solving 
the resulting differential equation yields 



l--i?(/,c) 

V 



l--/?(/o,c) 

^^0 



(A12) 



where 



R{f,c) 



(nX + 1 



2 2c 



f -2/ + ---(l + ( 

A n n \ A 



-nX 



(A13) 



{nX + l){nX + 2) p 

The equation ( [A12| ) is independent of the constant c by ( |A10b| ). We have introduced c as a 
trick for finding the final equation for all the trajectories. One way ( [A12| ) can be solved is 
to have 



l--i?(/,c) = 

V 



(A14) 



for all / and v on the trajectory. Different c correspond to different solutions and give all 
possible trajectories (either directly or as limiting cases). 

APPENDIX B: DETAILS OF TWO-LOOP POTENTIAL 






(b) 



(c) 




FIG. 7. Two-loop diagrams contributing to the effective potential. The heavy dots represent renormal- 
ization counterterms. Each line represents both heavy and light mode contributions. 

The two-loop diagrams (a) and (b) of fig. give the following contribution to the two- loop 
potential: 
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V, 



(a) 



^U{u + v)J{ml ml) + l-{n-l)UuJ{ml ml) 

+ ^{n-l)J\r{u + v)J{ml, ml) + ^{n-l){n-2)UuJ{ml, ml) 



V, 



where 



(b) 



' t I \2 2 2 2\ I ('^ -'^J 2 r/ 2 2 2n 



1 

12' 
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(Bla) 
(Bib) 



J{x,y) = J{x)J{y) , 

1 r d'^k 



and 



J(x) 



I{x,y,z) 



1 -p rf\ 



(B2) 
(B3) 

(B4) 



(27r)2'i 7 + x) (g2 + [(k + g)2 + 2] ' 
Contributions were considered from all mixtures of light (ma) and heavy {m^) lines, with 
masses given by (|3.5|) . The effect of diagrams (c) and (d) involving one- loop counter-terms 
is to replace J{x,y) and I{x,y,z) in ( pT]) by 



AfJ{x, y) AfJ{x, y) = AfJ{x, y) + - [xJ{x) + yJ{y)] , 
Af'lix, y, z) ^ N^i{x, y, z) = N^I{x, y,z)-- U[J{x) + J{y) + J{z)] . 



(B5) 
(B6) 



Including two-loop counter-terms corresponds to simply throwing away any remaining 1/e 
and l/e^ pieces in the potential.^ 

The general form for /(x, y, z) for arbitrary arguments in arbitrary dimension is given 
in ref. |T^. We shall need only the following special cases for the various expansions we 
make:0 



Note that the relationship between M and (/) has non-trivial dependence on e through M . Including the 
two-loop counterterms corresponds to throwing away remaining terms of the form A/'''^^e^"u'^m^*^^~'''0^° 

in V or e-«u3~a^2(2-a)^j2a -j^ J^^^v. 

Note that our e is 4— d where as that of ref. [|l5| is (4 — (i)/2. Also ||T^, there are some typographical 
errors in ref. jl^. In their equations (5.8-15), each explicit factor of J{w) in those equations (but not J(w, w) 
or J{v,w)) should be multiplied by 1/k. The factors of /i^' in (3.4-5) should be eliminated. The left-hand 
side of (4.8) should be x'^+y'^+z'^ . The second term on the right-hand side of (4.26) should be multiplied 
by 2. 
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(61 

N"^I{x, X, x) = x< — - H — (6 Ina; — 9) — 3 In^ X + 9 Inx 



-6^3 



6 



vr 

-In 2 

6 



(B8) 



where L{z) is Lobachevskiy's function, defined by 



L(z) = — / dx In cosx . 



(B9) 



and the value of interest is given in eq. ( |6.33| ). 



Ignoring the hght mass ma in (pl| ) and expanding to leading-order in e with u=—v gives: 



p^heavy 



(n-1) 
144 



J\fu^ (j? ml 



6 V/i2y 3 V/i2y 



(n+li 



+ 0{eV2] 



asym • 



(BIO) 



The coefficients of ((4 .81) and (|6.13|) and Cm2o of ( |6.9|) and (|6.10|) may then be extracted. 
If one makes the above expansion without assuming u+v is precisely zero (and expands to 
first order in = ^J\ffi'^(^u + v) (f)"^ + 0{m'^)), one obtains ( |6.17|) . 
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